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I. SUMMARY 


The two-scale composite rough surface model usually considered is one 
composed of large undulations over which small irregularities are superimposed. 

This general model may be further subdivided into two other models: (1) the large 
undulations are larger in dimension than that of the illuminated orea so that within 
the beam of illumination the picture is a tilted perturbed plane; and (2) the large 
undulations are of such a size that at least several undulations can be found within 
the beam. The first model is essentially the small perturbation model, since the 
effect of the tilt can be accounted for either by a change in the angle of incidence 
or by resolving the incident plane wave into horizontally and vertically polarized 
components, or by both. The second model is much more complicated and has been 
approached in most cases with a non-coherent assumption, i.e. the contribution from 
the small irregularities may be computed by summing powers from the large facets 
constituting the large undulations. The total contribution from the composite surface 
is then taken to be that from the large undulations plus that from the small irregu- 
larities averaged over the large undulations. 

If the non-coherent assumption is not made, the total scattered field from 
the illuminated area must be computed before evaluating the power which is the 
approach adopted in this study to calculate both the vertically and horizontally 
polarized scattering coefficients. The surface is assumed finitely conducting and 
homogeneous; the surface roughness may be non-uniform. To gain insight into the 
mechanisms of scatter, results are compared with those obtained by previous theories. 
The special form of the terms due to excluding the non-coherent assumption and the 
meaning of such terms are discussed. 

Based on Cox and Munk's [9] relotion between the meon square surface slope 
of the sea and wind speed, curves are drawn showing the backscattering coefficient, 
CTpp, as a function of wind velocity and of the angle of incidence. The comparison 
with NRL 110,11| data shows satisfactory angular agreement for both horizontal 
and vertical polarization, especially for incident angles larger than 30°. For 
smaller incident angles, NASA/MSC [121 data have been chosen for comparison and 
close agreement is again observed. 



II. INTRODUCTION 


Many scattering theories [1-4) using two-scale rough surface models have 
been developed in recent years. The general model is a large undulating surface 
with small irregularities superimposed. More specifically, two types of models 
have been considered: (1) the large undulations are larger in dimension than the 
Illuminated area and thus within the beam the picture is a tilted perturbed plane, 
and (2) the illuminated area contains at least several large undulations. For model (1) 
the problem remains essentially the same as o small perturbation problem since the 
effect of any tilt can be accounted for by a change in the angle of incidence and 
by resolving the incident plane wave into horizontally and vertically polarized 
components [3]. Such a simple treatment is not possible for model (2) for which two 
different approaches are in existence: (i) the non-coherent approach where a non- 
coherent assumption is used to simplify the problem [1],ond (li) the coherent approach 
where the said assumption is not made [4], The non-coherent assumption referred to 
here is the one defined by Semenov [1], i.e. the contribution from the small 
irregularities may be computed by summing powers from the large facets constituting 
the large undulations. This assumption implies that the calculation of the contribution 
from the small irregularities in model (2) is identical to solving the entire problem 
using model (1). Of course, for model (2), the contribution from the large undu- 
lations must also be computed and this is the major difference between the two models 
when the non-coherent assumption is made. Differences between models (1) and (2) 
are further magnified if the non-coherent assumption is not made. Thus, the total 
scattered field must be found by integrating the total field on surface before calculating 
the power. Hence, it follows that terms due to integrating the first order perturbed field 
over the large undulations within the illuminated area will show up and give an 
explicit indication of the interaction between the large and the small scatterers. 

Such an interaction is restricted to be an average operation in the non-coherent 
approach and this defines the major difference between approaches (i) and (ii). 

This paper discusses another coherent approach for model (2). The composite 
surface, If (x,y)= Z(x,y) +s(x,y), is assumed to be finitely conducting and homogeneous 
with Z(x,y) representing the large undulations and s(x,y) the small irregularities. 

Z(x,y) and s(x,y) are to be generated by independent, stationary, Gaussian random 
processes. 
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The approach is a modified Kirchhoff's method employing the equivalent 
surface field, i.e. the surface field on Z(x,y) estimated by the tangent plane method 
is modified to include the effect of s(x,y). Once this equivalent field is obtained 
the problem reduces to a single surface scattering problem, i.e. scattering from 
the surface Z(x,y). The concept of equivalent field was advanced earlier by Bass 
and Bocharov [5] for scattering from a single surface. Results obtained by this 
approach are simpler and reduce more readily to special cases of single surface 
scattering than Fung and Chan's approach [4], where fields on the composite sur- 
face f(x,y) were considered. 
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Ill, THE SCATTERED FIELD 


The far zone scattered field due to on incident plane wave on a rough surface 
Z(x,y) (Figure 1) is given by a special form of the Stratton-Chu integral (6], 



expC ng) ds 


( 1 ) 


where a time factor of the form, exp(jojt), is understood; is the position vector 
pointing from the origin of the coordinate system to a surface element dS;j ^2 is 
is a unit vector in the direction of observation; R is the distance from the origin to 
the field point; JH are the total electric and magnetic fields on the surface; 
k is the wave number in air; pis the intrinsic impedance in air, K^= -jk exp (~jkR)/ 
(4;rR), and ^ is the local normal to the surface. 

The basic problem for finding E is to determine n x E and n x H at any point 

* ah ah ah ah ' ' 

on Z(x,y). To do so It is necessary to set up a local coordinate system at the point 
In question. A possible set of local coordinates is 

1 = ( ~ + Z/ + I ) ^ (2a) 

5 = (ix-n,) / |ix<n,| (2b) 

1 = ^ 1 (2c) 

AH 


where Z^, Z^ are the partial derivatives of the surface Z(x,y) with respect ro x and 
y respectively, and j, k are the unit vectors of the (x,y,z) coordinates. From 
the definition of the local coordinates we see that the tangent plane at the point 
(x,y) on Z(x,y) coincides with the plane R-y. If the small irregularities were absent 
this would be an infinite flat plane. However, with the small irregularities present, th 
becomes a pei'turbed plane and the local scattered fields may be found by Rice's 
theory [7]. This has been done by both Rice [7] and Valenzuela [8]. Hence, let 
us assume that in the local frame 






t 



H = ^ H- 




iH- 

AH 4 


(3a) 

(3b) 
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2 2 

Assuming Z , 7 .^ \ , wc may rewrite the backscattered field in terms of 


the local field components as (see Appendix I) 


£s - " E^- tTiHa cos0)2j, 

/sLn 0'] + (icos 9 + k, sine) [( + -n Hy cos 6 ) 

-b ( Ey COS0 - 'a )Zy/sin e'j] expC-j]<t'X ) dx dy (4) 


where k, = kn,; n, = i sin G - k cos G; 2L, width of the illuminated area; G' is 
the local angle of incidence. ' 

For computing polarized scattering the Z - terms in (4) are unimportant 
and, therefore, reasonable accuracy may still be achieved by leaving out the Z^- 
terms. However, for depolarized scattering the field expressions are complicated 
by the presence of the Z^- terms and the local angle appearing in the denominator. 
In addition, depolarization due to the split of the incident polarization into locally 
horizontal and vertical components as a result of Z^ also complicates matters. In 
short, there is as yet no satisfactory method for estimating. the depolarized scattering 
coefficients from (4). 

For a horizontally polarized incident plane wave of the form 


E j_ = j exp [ - j k ( X sin 0 - z cos 0 ) J 

%V^ MM ' ^ 


there correspond two locally Incident waves (Appendix I). However, for polarized 
scattering there is no significant error if we take one of the local incident waves 
to be 

E = y exp [- j k ( ^ slue'- E cos 0' )] 

and ignore the other one. 

By applying Bass and Bocharov's [5] conoept to Valenzuela's results [8], 
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the equivalent fields on Z(x,y) up to the first order may be shown to be (see Appendix li) 


“ JJ i^vQ'EX dudv (5°) 

Ey = ( '+ R.) exp (-J Tx (u^bc)Q'EX dudv (5b) 

T\Hi = (i-Rx)cos0'exp(-jk,-r)+J| Tju\b-c)/k + kc | QEXdudv (5c) 
■i\Hy=// [uv(b-c> Q'/k] EX dudv ^ (5d) 

where EX = exp (-jux - jvy + jkZ cos G), 

Q ' = j ( k k^) S ( u - k si-nS’. v ) [ 2 -fr ( k^c + k'^b ) ] 

= I + Rx 

( k*- u* - V® , k^ ^ 

- j C + V* - k^ ) ^ < u* + ^ 

S(u,v) = two dimensional Fourier Transform of s(x,y), 

^ r Co, CCo) $ 0. Co = C k’^-u®-v^)'/* 

I - Co, J'm CCo) > 0 ^ 

and ym means "the imaginary part of." 

The limits of integration in (5) are from - “^to 

Similarly, for a vertically polarized incident plane wave of the form, 

E „ = ( ;^eos 9 + k sin 0 ) exp [~jk C-xsinO ~ z cos 0 )] 

the local fields on Z(x,y) up to the first order are 

Ei = C I - R„) cos6'exp (-J k,x) 

IJ T„Q'[bu si.Tv0'-( v2+bc) Ucos^'/k' J EX dudv 
Es “ //T/,Q'[bv sin 0’tu.vi< cosy '/ k‘] EX dudv 

=JJ T„Q'[uvic-b) cos ’p/k vk sin B' ] EX du dv 

= - Cl + R. ) exp(- j Ik, -j:; ) 

ku sin 0'+ Cv^c -/b-ck^)cos<A/k'] EX drxdv (6d) 

where T„ = 1 -i R„ . 


(6a) 

(6b) 

(6c) 
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IV. THE SCATTERING CROSS SECTIONS 


If we use a linear approx imafion for the Fresnel reflection coefficients 
(see Figures 2 and 3) and the local cos G' and sin G', then 

Rx = FC.10) + Rx’ Z* ; 

R„ == Zx : 

Rj! =i-£|< Rx( 0) stn 6 / ( k C05 ^ ; 

R == [ 2 k ( k' ^ ) Sen 9 ] [ k cos C k cos 0 + k cos ; 

COS Q' ~ cos 0 + Zx sm 9 
sin e‘ - Sind - Zx 9 . 

cos f = [ I ~ c k/k'f stn^B ; 


and (5) and (6) may be substituted into (4) to obtain the backscattered field (Appen- 
dix III). The scattering coefficient defined in terms of the scattered field by 

, \ 

Opp = ^-ttR ^j^sp'j^sp y / (2L) 


(7) 


can now be computed. < . . .> is the symbol for ensemble average; * is the symbol 

for complex conjugate. Some identities for ensemble average useful for evaluating 

are 

PP 


<Zx exp[ (z -z')]> = < Hx' j ''i cr-z')] > 

< Z, zl- expC jv, (Z.2')J>. 




-V 


G ( u' , v' 3 


as t 

d u. 




( V - V ) 


<S ( v') 5 ( u , V )* > = 2 TT cr,^ W ( u. V) <f ( u. - u.' ) (fi v- v'> 
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Comparison between llic reflection coefficient R(0, Z^) and its linear 
approximation for c>-/77r =5.61 (a) R,(0,Z^) 




Comparison between tbe reflection coefficient R(0, Z^) and its linear 
opproximation for €-/V/V=i.e\ <b) RiCe.Zx^ 
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I 1 I I \ I I I I I I 
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Z 

X 

Fig. 3. Comparison between the reflection coefficient R(G, Z^) and its linear 
approximation for C^= ZO (a) R„(0,Z^) 
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Fig. 3. Comparison between tlie reflection coefficient R(0/ and its lineai 
approximation for ZO (l>) Rj.(0,Zx) 
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where Z - Z(x,y); Z' - Z(x',y'); (x,y) and (x',y') represent in general two dif- 
ferent points within the illuminated area; (u,v) and (u',v') are similarly defined as 

(x,y) and (x',y') in the wave number space; rv = x - x', p = y “ y'; p = p(n'p) is the 

2 2 

correlation coefficient of Z(x,y); ‘-’j , a ore the variances of the surfaces s(x,y) 

and Z(x,y) respectively; 6 ( ) is the Dirac delto function and W (u,v) is the rough- 
ness spsetrum of s(x,y) related to its correlation coefficient p^ (o',P) by the Fourier 
transform. 

The general form of a may be written in terms of the sum of ai and 

pp • PP 

^^2pp (Appendix IV); 

(aL-loCIK2L-|t3|) f - 


ipp 


tppi 


— j ^ Vjj ( App Bpp A pp Bpp ^ 


aP 

aa 


- CT 


B 


Pp 


[ 


* ' t 

Bec" ^ ^ \ < act 




exp [ - j V, a - Vg ( I - f ) J d «. d ^ 


( 8 ) 


where 


V = 2k sin 0 and v = 2k cos 0; 
X z ' 


JsV fsi- 

2 J.pL 


■'zpp 


- 

''if 1^/2 

MtCppI 



^ d(x^ ^ ^ I ^/-ppi w 


ppl -||^f^e[TCpp(TDpp+R'Cppf J ) 


ex 


p[- j( u + k 5Ln0)a - j vp - K(l-j’)j da dvdos d(3 
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( 9 ) 


where Re means "the real part of" 

q = u - k sin © 

W = W(q , v) 

T = 1 + R 

P P 

Rp= Rj^ for horizontal polarization 
R,, for vertical polarization 

R ' = derivative of T with respect to Z 

P P X 

2 2 2 

K = 4k cr cos 0 . 

For horizontal polarization, 

hh " °lhh ""2hh (10) 

in which 

AhK = 2 COS e 

Bhk = £ CRx Sen 0 + RiCOSQ ) 

Ckk = Q(u^(fa-c)/k+kc + (u-^+bc.) cos d ] 

Dkh = Q C be) sLn 0 

Q = CR’®- k"') C 2-n t k'® b) 

For vertical polarization, 

^vv ~ ^ivv ^zvv (11) 

in which 

Avv = - 2 R„ COS 0 

= - 2 ( Rlcos0 + R„sLn0 ) 

Cvv = Q[[ciCb+ k cos fl ) j sCn 9 + [cos 0 C - v^b - c k^) 

- k C v^4 be ) j cos i> / k' j 

C>vv = Q |si.n0 ( ku sin 9 + ( v*c - v^b - ck^) cos ^/k' J 

+ k^sin0 cos 0 [cosS( v^c - v^b - ck^) - k ( v^4bc ) j/k' ’cos ^ 
- u cos 0 ( b + k cos 0 ) } 
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The expression for fT^pp identical to single surface scattering 
obtoined by the Kirchhoff's method. The form of cj^ may not appear 
because most of the cases discussed in the literature are special cases, 
the next section. The expression for cr2pp niore complicated; hence, 
to examine its meaning by considering special cases. 


results 
fami I iar 
as noted in 
it is best 
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where Jo ( ) Is the zero order Bessel function of the first kind. Eq. (13) is the 
bockscotter integral most often discussed in the literature. It is important to note 
that some of the conditions under which (8) reduces to (13) are loosely defined 
in terms of inequalities. Hence, the precise region of validity for (13) remains 
obscure . 

It is interesting to note that as L goes to infinity and for sufficiently small 
K (i.e. small a/x), (8) may be approximated, except for a specular-type term, by 

^1 pp ~ ^ 2 ^ I Appj J y ( ^ ) vJo ( 2 k ^ si.Ti 0 ) . 

Thus, 

e 

^1 KK “ 8 0"^ cos^0 Ri W ( £ k sin 9 ) 

^(vv = ® alecs'" q| WesksLuQ) 


16 



where '5 seen to be identical with the first order predictions of the small 

perturbation theory, |2,8| when the surface under consideration satisfies both 
assumptions of the Kirchhoff's theory and the small perturbation theory. However, 
this is not the cose with indicating that different approaches need not lead 

to the some results, because the degrees of approximation for different theories are 
different. 

If the small irregularities ore absent, oo =0 and a reduces to oi 

2pp pp Ipp 

as expected. 

If the large undulations are absent, i.e. Z(x,y) = 0, then becomes 


2L r2L 


iPP 


J 




k^lA 


2L -'-2L 
2 -2L 


(2Lr 


^PP 


' -JV, a 
2 da d f3 


PP I 


4 TT 


r ( £L-la| ) e 

“'-2L 




OC 


which is a specular-type term that behaves like sin x/x. 
With Z(x,y) = 0, CF2pp 


(2L-|ot|)C2L-|/3| 


spp 


^-2L-'-2L 


(ZL)‘ 


fce 

-LlJ 


TCppI 


W 


C ll - k si.n 0, v) exp j ( a + k sen du dv doc d / 


If 2L can be taken to be infinity, CT 2 pp reduces to 


Thus, 


2 r“ r® s 

Tjpp = 2('7rk0"^) / I |TCpp| W(w.-kSLTi0, v)i(u+kSL7i0)d’cv) dud 

J-CO J- 03 , , 

(14) 


vv = ® ^ I R//COS^0 + t/( k' - k*) si.n*0 /2 k' ^ I 
'Weak sLn 0,0) 

^ 2 kh = 8 k"* cos^O I Ri.1^ W (2 k si.u 6,0 ) . 


(15a) 

(15b) 
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The above results are identical v/ifh the first order results obtained by the small 
perturbation method |2,8| ^Appendix V). Note that 2L must be taken to be infinity 
for 2pp reduce to the perturbation results because the mean plane for the 
perturbation model is an infinite flat plane. 

Let us now examine the first term in <n when Z(x,y) is not zero and v/hen 


2L can be taken to be infinity. It has the form 


'2pp 


^ ao c<o fco , 2 

X,X.LLlTCpp| W(u-ksm0,., 

exp j Cu-t k si.n. 0)<x-jvp-K(l-f)3cJudvdaj.dp ( 16 ) 


The variables u, v may be Interpreted as the wave numbers or frequency components 
for the surface s(x,y). Comparison between (14) and (16) shows that in (14) only 
a specific pair of u,v values is allowed whereas in (16) all values of u,v are required. 
This means that the large undulations are responsible for making all frequency 
components of s(x,y) effective In the scattering process. They also define the 
appropriate weighting (through K and o(n,p) ) on the. contributions of the different 
frequency components of s(x,y). Similar statements can be made about other terms 
in <^2pp fXiey vanish with Z(x,y). Thus, cr2pp define expli- 

citly the interaction between the large and the small scatterers. This interaction 
vanishes when either s(x,y) or Z(x,y) Is zero. Scattering theories with the non- 
coherent assumption have this interaction replaced by averaging (15) using the slope 
distribution of Z(x,y). (The dependence of (15) upon the slopes of Z(x,y) arises 
when the incident angle is taken to be the local incident angle.) Hence, in such 
theories the nature of the interaction is assumed rather than calculated. 


If the 0,3 -integrals in 


converge fast enough so that within the region 


^_2pp 

of convergence 2L» o,p and if, in addition, edge effects are negligible, then 

■ 2L rSL roo r«> 


a 

app 


I 2 .^2 


/ CL pCL rco r 2 

/ / / if iTCpp -t G CTDpp R'Cpp ) I w 

-2L -'-2L -'-OO J-t» 


-t V, 


Re(GTCpp[TCpp GCJDpp- R'Cpp)]*) 

r'/> ) 


~T ^ C w ^ ^ ^ 

T G- C pp — 


exp^- j (a+ k sLnO )a-jv(3-K(i-j>)j dudv da dp 
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(17) 



where G = (ui k sin G) / v^. 

For isotropically rough surfaces, (17) may be further reduced to 

||TCpp-^G(T Dpjj + R Cpp ) I W 

Re (TGCpp C TCpp ^ GCTDpp + RCpp) ] *) 

+ 1 TGCppf (^y(u+ksLnej"4 

• du dv ? . (18) 


To compute the backscattering characteristics with (13) and (18), p(?) and 
W(K) need be specified. As an illustration, we assume that p(g) is Gaussian and can 
be approximated by the first two terms of its series expansion about 1=0 and 


W ( K ) = ( i V2 ) [ - ( K i 2 f ) 


(19) 


where J?is the correlation length of the surface, s(x,y). Under these assumptions 
(13) and (18) become (see Appendix IV) 

<r,pp = (8 m*cos®0 |App + Bpp taiaO] exp [- tan*0/( 2m® ) j 


(20) 


^2pp ~ ^ ^ 


:'0I{UCpp^GCTD„.R'c„) 


w 


+ Re([TCpp + GCTDpp+R'Cpp3] T*Cpp"G ) V, j 
Vj m ^exp|^-[(iA + ksi.n e/+ 2 ) j duL dv^ (21) 

Figures 4 and 5 show the general angular behavior of and for 
different values of the rms slopes of Z(x,y) and kpj of s(x,y). Since the major dif- 
ference between this theory and other scattering theories lies in 02pp' ^'9’-"'®^ ^ > 

‘^2hh ' plotted using (21). First order results from the small perturbation 
theory for a single surface given by (15) are also shown to provide a basis for compar- 
ison. In Appendix VI, all the identities used are rewritten for ease of reference. 
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in 



Q in degrees 

Fig. 5. Angular behavior of <T|^|^ as a function of the rms slope, m, and kcTj. 
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VI, COMPARISON WITH SEA CLUTTER MEASUREMENTS 


As an illustration, only tho results corresponding to an isotropic rough surface 
with Gaussian surface height distribution will be compared with the experimental 
sea data. The assumed surface model although not realistic for the ocean surfoce, 
is oble to predict the correct angular trend of the bockscottering coefficients for 
both horizontal and vertical polarizations with surface parameters of very reason- 
able sizes. 


An examination of the bockscottering coefficients, given by (20) and (21), 
indicates that the rms surface slope m of the large undulations not only affects the 
Cipp term, but also the cr2pp clear that the large undulations 

influence the returns from the small irregularities as mentioned in the previous 
section. To suppress os much of the effect from shorter waves as possible, Cox and 
Munk's [9] mean-squared slope measurement in the presence of oil slicks is adopted 
to estimate m for different wind conditions (see Appendix VII), 

It is generally agreed thot the surface spectrum of the small irregularities 
should vary like BK ^ (K = surface wave number and B = constant). In comparison 
with experimental data, the value of kl of (19) is assigned to be 2 (Figure 8) so 
that the correct angular behavior of the Gaussian spectrum approximates BK 
well over the angular range, 30° 0 70°, i.e. BK ^ is approximated by (19) 
with K = 2k sin G, the Bragg scatter condition. To bring the level into agreement at 
0 = 60°, we multiply the Gaussian approximation by the factor of 35.3 . Since a 
complete information of the increase in intensity of the high frequency part of the sea 
spectrum is not yet available, the wind dependence of cannot be uniquely deter- 
mined. Oceanographic investigations indicate that the values of B lie in the interval 
4.6 X 10 B5 l3,26 x 10 ^ [9,13, 14]. This implies that ka^ should lie in the range 

from 0.067 to 0.2 when BK ^ is equated to (19) at 60°. These values of ka^ are 
consistent with the assumptions of the small perturbation theory. 

According to the above arguments, comparisons of computed [(20) and (21)) 
and measured (NRL) backscatter characteristics of X and L bands for wind speeds of 
11-27 knots are shown in Figures 9 through 16. It shows fairly good agreement, 
especially for -incident angles larger than 30°. Since questions have been raised about 
the accuracy of the absolute levels of the measured scattering coefficient curves from 
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Fig. 11. Comparison of computed and measured backscatter characteristics. 
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SCATTERING COEFFICIENT (dB) 



Fig. 12. Comparison of computed and measured backscatter characteristics. 
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SCATTERING COEFFICIENT (dB) 



Fig. 13. Comparison of compufed and measured backscatfer characterisfics. 




SCATTERING COEFFICIENT (dB) 



Fig. 15, Comparison of computed and measured backscatter characteristics. 
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different missions [15], the main intention of these comparisons is to show the angular 
agreement rather than the agreement in absolute values. However, to ollov/ reference 
bock to the measured data, all levels of the data have been raised by 6 db rather 
than arbitrarily adjusted for each wind speed to obtain better fit to the theoreticol 
curves. Values of koj and m used in the calculation are as follows: 


koj 

m 

corresponding wind speed (knots) 

0.1 

0.09 

11 

0.13 

0.1 

14^16 

0.17 

0.12 

23-27 


It is noted that the value of ka^ increases with increasing wind speed. This obser- 
vation is in agreement with recent studies of the sea spectrum [16,17]. For smaller 
Incident angles more data points ore needed to define the angular shape of the 
CTpp curves. For this reason NASA/MSC data ore chosen in Figure 17. Agreement 
is observed between measurement and theory with m determined by Cox and Munk's 
clean sea measurement. More experimental data are needed to explain the discre- 
pancy in using slick and clean sea measurements. 

On the basis of the above results, it appears that it is possible to determine 
the wind dependence of the scctterometric parameters m ond ct-| • With sufficient 
experimental data these parameters may be determined more precisely for different 
wind speeds. 
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VII. CONCLUSIONS 


The results of the present theory Indicate an explicit interaction betv/een 
the large undulations and the frequency components of the small irregularities. As 
compared with ^^]pp/ *^^2pp represents this interaction decreases more slov/ly 

with the increase of the incident angle. This offers another possible explanation 
for what has been called the diffuse scattering portion of the angular curve.;-. 


38 



APPENDIX I 


THE BASIC SCATTERED FIELD EXPRESSION 


In this appendix, an expression of the backscattered field in terms of the 
local field components on surface is derived. The star'ing point is the modified 
Stratton-Chu integral: 

lie X X E - TT H)] exp ( j l< t;- ne )ds. (M) 

Assume 

E. = S y ^ zEi 

J *■ 

H = i > 3 Hy ^ iHt 

where 

X = / l2i I 

1 = ( I :ri. )/Do 

« i. 2y COS 0 + J (sin 9 - Z^cos 0 ) + J< sLnO J £)j' 

■c It z/ 

-? = y X E 

~ ( X - i COS 9/D. + Z, )( ( -V z! - Zy 

"lx ^21' I ~ [ ^if '*' ( “ 2 * cos 0 / j ^ 

2 2 

If «1, the local unit coordinate vectors may be approximated as follows: 

X “ X ' ^ CO5 0 / si.n 0 + is. Z, 

X ~ X ^ a ^ i ii sLn 9/sLnd 

L ~ 0-2) 
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with 


sin 0 ~ si-n 9 - cos 0 


For bockscottering (i.e, = H^sinO + i^cosG) tho following relations may 

be obtained 

~ 1 " i 

li.' * 5. “ ~ ( i. + Jt sinQ ) Zy COS fi /sLnG’ - j cos 0 

21’ * 1 = ^ i. ® Ji si"^ 0 ) ■ j 2y / slri 9 ’ 


with 


cos 0 « cose + 2 . sen 9 




■ ~ < i cos 0 t j< sin 8 )(£«■*■ Zy cos 9 Ey /sin 6 ) 

4jtcos9 Eg - HyE^/ sl-n 0 ) ] , 


(1-3) 


Similarly, 

2ia*C2i2’‘<21^if>}== ( ^ cos 6 ■r k. si-n 8 ) c cos 0' Hg - 2 g H^;/sin 9 } 

- j ( -t Z u cos 9 Hg / sin 0') . 

Notethot 

( i sine -J<C05 ) ^ C c cos 9 + k sin 8 ) = - j 
n I « i = ( -i cos 9 + k sin 8 ) 

A * 

Substituting (!~3) and (1"4) into {1~1) we get the bockscattered field as 

£ s = ~ K„ I i ^5 0 - ) - ( E,' t cos e ) Zy/slna ] 

+ COS 0 + J< sin6)[c E;x + ■>jHgccs0')+(Egcos0-'>]Hj)Zy/sine ] 

■ exp (- j Jk, • jr ) d^f di/ 


where 


(1-5) 


k, =i k n, 
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By using (1“2) it is possible to find the local incident fields. Thus, for a 
horizontally polarized incident field of the form j exp(-jki ' r), there correspond two 

,4^ /vU. ^ 

local fields, i .e . 

... , - Zy 

j exp ^os 0 -r i. sen 0 ) ^ir\Q‘ 


■e>‘p jii' *z: ) . 


(1-6) 


Similarly, a vertically polarized incident field of the form (^cosO + J<^sinG) exp 
(-j k.| • jr^ may be decomposed into two local fields, i.e. 

( COS 0 + Js sLn 0 ) exp ( - j is .1 ■ D ) 

~ ( -X cos 6 + i sbn 0 

(1 

+ 5 ( /sin 0 ' ) exp ( - j k, • r ) . 
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APPENDIX II 


FIELDS ON SURFACE Z(x,y) 


A. Horizontally Polarized Cose 


From the small perturbation theory, the total scattered E-field up to the 
first order in space may be written 

) ~ ff "Fj. av Q Ex dadv 

^ “X I 

E5 t ) =[expcj k Z COS0') + Rj. exp ( - j h ^ a')J 

■exp (- j R X sinB' )-fJ' Tx ( cb ) Q' E du. d' 


(ll-la) 


= E^o +Ey( 


(ll-lb) 


where 


= exp (- j u X - j V y - j kcosfi ^ ) 

Q' = j ( k'^- k*) S C «• - k sEri 0', V ) C 2 TT D ) 

D = U^c -4- k'^b 
Tx = I + R X . 

To obtain the zero order E-field on surface Z(x,y) we may apply the tangent 
plane approximation. Thus, in local coordinates its value at any point on surface is 

EyoCO,0,0) = I+Rx 

To obtain the correct phase relationship between points in(x,y,z) coordinates, the 
total zero order E-field at a point (x,y) on surface is expressed as 


^ yo ( , y j Z ) = ( 1-t R^) exp ( - j ill ' ) 


( 11 - 2 ) 
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The form of the first order field is not similar to the zero order in that it does 
not contain the incident field. Thus, taking as an example, we find the value of 
the field on surface in local coordinates to be 

E ^ , Co,o^o) = J j Tj_ U.V Q du d V 

Following Bass and Bacharov, the total first order field on surface in (x,y,z) coordinates 
is 

/■ r , .;ux-jvy-jkzcos6 

du.dv_ (II” 2) 


Thus, the complete set of fields on surface in (x,y,z) coordinates is 


Tx uv Q dudv 

(Il-4a) 

Rx)expC-iiSi - E ) 

(ll-4b) 

-// Tx (^*+bc)Q'EXdu.dv 

9'< (- Rx ) exp c-j ) 

(ll-4c) 

Tx [ u (b-c)/k + kc ] Q^EX dudv 

rj^fuv<b-C>Q^/k j Q EX dudv 

(ll-4d) 

- j V y - j k £ cos 9 ) . 

(11-5) 


where 

EX = exp ( - 3 

B . Vertically Polorized Cose 

For a vertically polarized incident plane wove of the form 

E = C'i.cos0-*jksi.T\0)exp[^-3l'C'^SLn0-2 cos 0 ) j 

the local fields on Z(x,y) up to the first order are 


( 11 - 6 ) 


E 5 = cos 9 ( I - R/, ) exp ( - j ii' ■£ ) (ll”7a) 

+ 1^ T„ q' [ bu SL-n 0'- ( v®+ be) U cos i>'/ h ] EX dixdv 
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where 


~fj t„Q v(u.Cc - b)c,os^/i,’ - U sine'] EX dud ^ 

= - c ' + R„) exp ( - j k.-r ^ 

//^"^ (iruslneV (v®c -v^t-cx*;cos}iS//{'] 
du. dv 

E,, = exp<- j U.X - j vi| - j kacosQ ) 

Q = jC k - k )[a7r(k'% t k^c)] 5 C u.- ksLn 9 , V 
X, ^ - R,, . 


^ai-7b) 
(II- 7c) 

(ll-7d) 
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APPENDIX III 


INTEGRANDS FOR THE POLARIZED FIELDS 

In order to moke use of the formulas for computing the scattering coefficients 
shown in the next appendix, it is necessary to get the integrands of the field 
expressions in the right format. 

The general field expression is of the form 

E = K JJ( INT) exp (- j dy (Ml-i) 

where INT for the two polarizations are given below. 

A. Integrand for (INT)j^l^ 

(INT)^^ = COSe'Ej - 

= CO50'([1+ Ri(e’)}exp(-j j<,-r )-/jT,(0')(u4bc)QEX dadvj 
~casd'[l - R^(6')]e?ipC-jhrj:)-JfTi(d')[u-'(.l>-<^V^ + <^c]Q'BX dudv 
= 2 [r.( 0) + Ri'Z,] ccose + Z,sLn0) exp (- j Ji,' -j: ) 

- // TAe)[u.'cb-C)/k + l<C t Cu®+bc.)cos0 + 2^sLrteCu.^bc)] 

•Q5<0’)EX dadv 

— 2 ^ Ri(0 ) COS0 + CRi( 0 )stn0 + RjlCOS0] Zx] 

exp ( - i ) - J/ R1 Zx 1 

r ^ ^ ^ , (111-2) 

[d.kkT Djkk ] [5C9) + S Z, )EX dadv 
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where 

3(0')== 3 C0 ) +3 ' 2 ^ , TxCe') == Tx (0) -*■ Rx 

Ex = exp ( - j - J ''a -r j h 2 cos 0 ) 

3 = = C ^-S/ac u. - (< Sun 0 ) 3 h cos 0 

D|Kk = [ a ( b- c)/l< + kc t (u*+bc) cos 0 ) Q 

DaKK =(a*-fbc) Q si.n.6 , (2 = Q-5(0'3 = QS (w- ksLn0> ) 

Q = j ( k' k^ ) ( 2 Tr C U c + k'"b ) ] ■' . 

Linear approximations hove been used to rewrite results in terms of the incident angle 
rather than the local angle. 

B. Integrand for Evv 

( INT Hy COS 0 ' 

= COS0'[|- R„(0')]exp(-ii5,i'j:) + J/T,,ce') Q' 

[ bu. sin 0 - C v®+ be ) k COS } EX Ju. dv 

+ cos0'{-c 1+ f?„ce')]e^pc- j b Z ) 

+ q' [ kusin G'+ (V^C- v^b - c.k’^) cos /k ] 

■ E X du dv J- 

= - 2[r„( 0) C.OS 0 -r I; R,' cos 0 + R„(0) Si.n0 J Z, J 

• Ji' *Z ) + // T. (9 ') QS(0') 

l^u. ( b + k COS 0 ) sin 0 '+ [cos 0 ( v*c- v*t - ck*) 

- < v*-» be > k } cos EXduudv 

+ J J T,^ (0') (J 5 ( 6') sin 0 [ ku.sin0 + ( v^c- - c k*) 

■ ’f‘/k' ] EX du. dv 

= - 2 I R,AQ) COS0 + [ r' cos 6+ R„(0) SLT10 J 2^ I 

• J Jl'-Z) -t//[:T,.(0) + RjZ, J[D.vv-^D,.v2, 

[S(G)-tSZ,j EXduudv 
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(II1-3) 





D WV 

DzW 




( 6 -.- 


9 ^ 


e> 






,/(k''co^+^ 


^^005 


e 


+ e» 

^coS. 


e ) 


. = U^>c.o& 

Qcos« . e^ C'^<^ 

&' , etk'^^'-’^ 

. Q &tn ® I 

C, ' ^ 


Cv 


4 ^/ 


<')• 
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APPENDIX IV 


EXPRESSIONS FOR THE SCATTERING COEFFICIENTS 


A. Consider a Field Expression of the form 

Ei=Ko|/ C App + Bpp Z, ) e d-x du 

-L ^ ' (lV-1) 

the product of E, and its conjugate is 

E, •E,“‘ = | Ko|* //// ( I App|^+ App Bpp Z«'+ I Bppl^ Z, Z«' ■+ App Bpp 

-L 

exp C -j V, cx-x') + j Vj5<z- z')]dx dy do<'dy', (|v_2) 

Taking ensemble average of the product yields 

<E, -E, >= /Kol iLlLlLo'-f-i-x' ( ~ ^ 

• - 1 Bpp l" ( -ip + j} e^p[-i v,oc- v|o-*0-f )] 

da d^ d X ' d jj ' 

= l(^ol J J ( SL " I I ) ( 2 L - 1 p I ) I AppI - J ( App Bpp ■» App Bpp ) 

(IV-3) 

where a = x“x'; /3 = y"y'; A^^ and are not functions of «, ^ ; K = V 2 ffj af=<HZ>, 

The scattering coefficient ^.jpp corresponding to the field given in (IV-1) is 

^ipp = ^Tf R*< E, •£, > /(2L ) 


^ x:: i:: {■ a„i'- j ca„ , a,*, b,p > 


•exp[- j v,a. - K ( I - f )] 


(IV-4) 
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If the od , f3 -integrals in 
convergence 2L»«, ^ , and if edge effects are negligible, then 


O’, converge fast enough so that within the region of 

IPP 


_ k' 


ipp 


-4 Tr 


^2Lr2L - 2 w 


pp App Bpp 




->pp 


•expC-j'^K*~ B (I - (3 

= i |App+-r Bpp|%xp(-J''-^-K('-?0 d/3 . 

(lV-5) 


For isotropically rough surface, (IV~5) reduces to 


'PP 


r 2 *- 

I ( ^PP 

^ /•% 


B 


pp I 


Jo ( ? ) e 


KO-?) 


t 


(IV-6) 


where Jo( ) Is the zero order Bessel function. 

The identities useful for getting (IV-5) from (IV-4) are 


(1) “S' expC-j v,a- || exp [- j V,oc - K(i-f )} 

+ - B(-H)''}expr-jV,oL-KO-f)} dot d /3 


2 L 


d/3 


.^L 


(2) 


J J ^ e\p[- j v^oi- K(i-P)jda.djQ= J^'~ ~ exp[- j V,a- K - ?)j 


2L 


'-2L 


JL^ 

K 


f exp [-jV^OC- doi d|3 

-^- 2 L J-ZL 


B. Consider a field expression of the form 

Ez = II ( II < R 'H^)(D.+ DaZ,.) (S-S'Zx) EXpdudv Jdxda^l^_^^ 

where and D2 are both functions of u and v; S' is the derivative of S with respect 
to Zx; 

EX p = exp [_~ j ( u. + k sLh 0) X -j''y + j 2 k cos 0 3 . 
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It is possible to write Eo same form os (IV- 1) and identify A and B v/ith 

2 ' ^ PP PP 

Zx -term being ignored. EXP may be taken to play exactly the same role as exp 

[-jv^x -I- jv^z] without affecting the final result. Note that appropriate subscripts 

should be attached to D^, D 2 , T, and R, depending upon the polarization states. 

However, these subscripts have been left out here for simplicity of writing. The 

basic form of the scattering coefficient for this field is again given by (IV-4). The 

corresponding coefficient terms can be shown to be (see section C) 


< I Appl"> = 2Tr 0-' // iTD.I^Wdadv 

- 00 


(IV-8a) 


< App Bpp > = 2 - T D. C ( T D. R'd, ) - W -t C t D. ■^yv-g'br 

C I Bpp|^> = fl {iTDa-^ RD.l^W + ) TD, 

Re [tD.CtDs 1 - R'D, j I da dv (IV“8c) 


+ V, 


-CC> 

3W 


where < . . .)• is the symbol for ensemble average performed on s(x,y); Re means 
"the real part of; * is the complex conjugate sign; 

W = W(q,v) 

q = u-k sin © , = 2k cos 0 , 

Although (IV- 1) may be used for (IV-7) , it is more convenient to write the 

scattering coefficient for the field given by (IV-7) directly in terms of and D 2 

instead of A and B . Thus, 

PP PP 


zpp 


= • 4 .TT R /(ZLT 

- J 7 ( ij 


k <r, 

3 


-EL 


TD.ppI W 


V, ^ Re (TD.pp[(TD,pp-^ RD.ppf w-v ^ T D,pp^'-|^ ) 

V 


(iTD.pp-t R'D,pp|"w- jTD,ppl^^'^ 


R«^T'D.pp(TD^pp + R'D.pp/j) (iv-9) 
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If the a, fi -integrals converge fast enough so that within the region of convergence 

2L>>«, ^ and if edge effects ere negligible, then for isotropically rough surface, similar 

to O’ , , (IV-9) reduces to 

Ipp' 


^zpp ^ fo J-a> I ^ ^ R Dipp ) 1 W 

+ fy Re (QTD.ppfTD.ppt G (TD 3 pp-t RD,pp) j""} 


"1^ iTGD.ppf j y C ^ Jc ut kSLne + v" )dadv 
[ - K(l-f ) j ^ d 5 


(IV-10) 


where G = (u + k sin © )/ v^ . 

C , Consider o field expression of the form 


E = ( AoS + Bo S ) E, EX P dudv dy , 

-L. *oo 


(IV-ll) 


then 


E •£*= |K,|' ;/// I jpf [ /AJSS.% A.6.*ss;'-. A.*B.3*S: 

^ I S,' ”] £, Hv 1 EXPI* dll dvd Adv' j d* /[y_| 2 ) 


where Ao, Bo are both functions of u, v; S - S(u,v) and S-j = Si(u',v'); EXP = 
exp [“j(u + k sin 0) x-jvy + jk cos 0Z). The following identities of ensemble 
overage over s(x,y) ore needed for getting (IV- 13) below: 

<SS*> = 2'rrO’,^W(u-U sL-n. 6,v) SCa-a.') ^Cv - v') 

< S s!*^> = 2 TT cr,^ — k cos © S(.u-- ll' ) S k V- v' ) 

< -S S. > = 2-TT k^ COS^0 ^ iiVL-iX ) ^ L V- v' ) 
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to 

■ ~^L 

""’"withr. 

^c^.'6 "’'^f//nA.,v 

<■ + /»„ e^ -j ,.av , -oo ' V 

•"t -s^r-. K( 4 if, , f 

<,ft a, ^ ”'■"’ "V- 3 J end (,v-,, 

'®pp'> = £>0-' f7r,„ a »e (,000 

J J t IA^I Vv/ + *rj 

■"'e.l%^os'e ^ j * -^T^kooos 

=0 ^ e e J 0 / V 

j ofuc^v 

*" "■ «V- 8 c) »,•„, 

60 ::-- OjT, ^ 


52 



If, in addition, we assume that autocorrelation of the large undulolions to be 
Gaussian and integrate out the Bessel function of (IV-6) and (IV-10) , the following 
expressions result: 

‘^ipp= T (4 no.^Cos'^e)'' ( App-r BppWe 

OCr ^ 

^epp~'nl< cT'i ||Cpp + QDppj \a/ + RetcCpp~t-Q- Dpp ) Cpp G- J 

-oo 

(lV-16) 

0 9 

where W^)= (j? /2)exp [-(Ki’/2) ] is the roughness spectrum of s(x,y) related 
to its. correlation coefficient by the Bessel transform; Jl is the correlation length of the 
surface, s(x,y); f (o) is the second derivative of p evaluated at zero; 

'»n*= cr* I ?"(o) 1 

W = W ( v^ ) 

Cpp = T D ,pp 

^pp = ~l*Dapp + R Dipp , 

To obtain (IV-15) and (IV-16), we have used the approximation and the 
identity os shown below 

0) exp C- K C 1 - y ) ] ~ exp [- K I ?"Co) | / z ] 


( 2 ) 


X J. < V, ^ ) exp ( - a ^ ^ 


!^Gx p ( - a. ) 
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we get 


Comparing {iV-]), (!V- 7 ) with (il(- 2 ) ond (l!!- 5 ) and using (iV- 4 ) oncl (iV- 


9 ) 


0 ) 

with 


( 2 ) 

with 




AtiK “ £ R JL ( 6 ) cos 0 

Bkk = 2 fRx(e) siu e+ cos 9 3 

f<c+(u^+hc;co59] 
Dzhk = Q (u^+ be) stn 6 . 


vv 


= tr.. 


cr. 


aw 


Bvv" ~£[R,, (.©) slu 9 + R„ cos 0 ] 

Divv ~ Q|[u_cbtk cos 0 ) ] slrv0 -t [cos $ Cv^c 
-v^i> - c - ( v^+bc) k ] cos ^0_/|< ' I 

Dsvv ~ Q ^ slug [ ku.si.n0 f Cv'c-v®i>-ck^) cos ) 

+ k^sc-n-e cos 0 [cos 6 ( v^c- v®L - ck^)~ hc)k 

/ ( k' Vos 5>5> 3 - u cos 0 ( b + k cos e ) j 
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APPENDIX V 


IDENTIFICATION OF TFIE DIFFERENT FORMS OF THE 
SCATTERING COEFFICIENT FORMULA 

To compare in (14) with the corresponding Valenzuela's and Wright' 
scattering coefficients, we rewrite (14) in the following form 


= 2 ( 


2 

! ) 1 W(u-ksln6,v) 

CO 


i(a + k5i.n8) ^(v) du. dv 
2 kcos 9 I 

" k'cos6+ kcos«)> I 2L I * 

Cvv = Q(a(bf Ucos6)sun9 +(cos0(v®c -v b-ck ) 

- k ( v*+ be) } cos “k/k J 

Q = ( k*- k*) (air C k*c + k’*b ) ) 

cos^<)> = I- k* si-Ti^e/k’* . 

By the property of the Dirac delta function it follows that 
u = -ksinG 


T./Cw(-ksin0,o) = T„Q[-k sLtiG ( kcos0 + kcos0) sin 0 + 

[cos 0 ( - k'k*cos 4 > ) - k cos 0 cos4> 3 'k/ k j 

= T,^Q 2 k^cosG (si.n*0 + cos^<k ) J 

= r c _ os _ e j f ] f- 2 k"co 50 esm ^0 + 1 - 

kkcos6 + kcos4> ^2irkk (kcos4> + kcos0 ) ^ ^ ' 

= -.2- Ck'^-k^)kcos^e(k'V(k'^-k^)5Ln^0 ) 

k^(k’cos0t kcos 4" )* • (V“3) 

Substituting (V-3) into (V-2) yields ^ 

^ - 8kV'cos‘0 ^(^ks.ne.o) 


2 

0-3^^ =2(^k<r,f |t„CvvC -ksl7x0,O)| WC2k5Ln0,o) 


(V-2) 


O' = 8 k cv; cos 0 

2W 


k ’’ * ( k cos 0 t k cos >p ) 
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(T given by Valenzuela and in his notation Is ( ~ 8 TT cos 6 ^ Psw^ 


( ^vv }val. 

S nr CO5 0 < Psvv> 

2 


4 'u (3^ COS^9 

(e - 1 ) [^ ( 5<-n^0 1 0 “ S(.n^9 j 



( e.cos6 + «/t- sLu^e )■' 

(V-5) 


, W ( £ k stn 0^ < 

^ ) 


The connections between his notations and the one used in this paper are 


/S - k 



(V-6a) 

e = k_;Xk^ 



(V-6b) 

2 

2 . 

W 


(V-6c) 

W “ TT 



Substituting (V-6) into (V-5) yields 


( <T 1 , = 4 -TT k C05 0 

' VV^VO.1. 


,2 1 V. 

( + ' ) - s-v-n^e J 


( Y- 0050+/-^'- sLTi^e 


2 CT, 


K' ( 2 k SLTl 0,0) 


= 8 k O', cos 9 


( k' ) C k'V ( k^ ) S>lyi^Q J 

k'^C k COS0 + k cos<^> ) 


(V-7) 


W ( £ k si.n 6 , O ) 
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whicti is identical to (V~4). 

Using the notations in this paper, we con rewrite Wright's scattering 
coefficient as follows 

(‘^vv)w 2 

4_2,„ „.^ee ^T:Cl<'-k")stn^0/(2l<'')| W (2 U s.’-t G^o) 


where 


= 8 k'‘cT,*| R..COS' 

_ k'cos9 - Ucos<^ 

" k'cos6 + 

Substituting R„and T„ into g^^yields 

g = R„cos^0 


(V-8) 


_ (kcosQ- kcos^) cos ^6 ^ 4 |<'®cos^0 ( h'~ k*) sLn 9 


k'cosG+ kcos4> ■ 2 k'*(k'cos0-t k cos<#>/ 

- _(k'’‘cos^e-k^cosV)cos^8f 2cos ^6(k'^- k^)sLn^e 
C k'cos S + k cos )® 

“ ('k'cosellL,i,)‘ f '' ® siu'sck'*- k’)] 

' cEoksTk = kVn'8 ] 


COS 9 


( k'cos 0+ kcos<#> )* 


[(k + k si.71^0- k'^sln^S) — ^(k’ - k^scn^G + k sLn^6)j 


cos^g Ck'^- k^)(k 4- k sln^G - k SLn^e ) 
k'*( k'cos 0+ kcos4* 


Thus, 


( 0-vv ), = 8 k^a."co5^0 I 

I k'\k'cos0t kcos4>)^ 


which is again identical to (V-4) . 

In conclusion, when Z(x,y) — 0, we get 


(V-9! 

2 

WC2 k 517X0,0) 

(V-10) 


cr. 


2VV 


^ ^vv)val ~ ^ *^vv^W 


= 8 k^o-.^j R„ cos^0 -t T,,^ k'^- k") 


3, , ,S , Z, SC7X 0 


ak'^^^l W( 2 ksf.n 0 , o) . 

(V- 11) 
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APPENDIX VI 


COLLECTION OF IDENTITIES 

All the identities used in this report are rewritten in this appendix for 
ease of reference. The order of appearance of these identities does not correspond 
to that in the report. All the integral identities are given first and then the iden- 
tities for ensemble averages. 

I f ^ f 

s ( u., v) = 7;;r J ^ SC9<, y ) exp(- jTox-jjv) dx 

2. J| f, C dri 

3. If W (|S3 =|"X( ^ 

then W'(/33= =-|”j. C/3 O 

a?, 

Jo 3 ^ J|C/3^)expCb^, :i:-|3Wc/5) + edge effect term 


^“^'expe-ax) J«(bx)dx 

= 2 a(ab)“r (a^-4-)/|; y:c (a^ 

Re <x>— 1^ Re<^^l Im b 1 ^ 
where P Cti+i) = -nPL'n^) if TL > o 

r (-i) = 

5. r exp ( - J« ( b^c ) dx 

= [ bV(2af*' ] exp C- bY^a,) 

^Re ^ a>o , Re cx>-l^ 
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6 . 


7. 

8. 
9. 


exp Cj hcas-x) cosn^K ci-x = j^"frJ^Cb) 

[ SLTi 0 exp (± J9< si.n a) cje = ± j S'" 3. CTC^ 

•'o 

exp(j'’‘sLn0 + jTnQ)d6 = (“’) C'X ) 

exp (tj'X cos Q t j ^ sin 0 ) d0 = 2 'rr Jo C Jx * ) 


10. J exp (jWt)doJ - 2 'JT Set) 

11 . < zz'> = 

12. < s s'> = (T,= ?, 

,3 < 5cu:v'y Sca.v)"" > = 2T- <5-.* W(u,v) S(U-u') SLV-v') 

<Zx exp[ jV^(Z-2')] > = < zl’ expCj ViCZ-z')] > 

= - i ^2 It ' - ? ) } 


where oc = -x - •x ’ 


15, <H,Zv «P[jVa(Z-Z-)J>..o'(;0 + <T=v/(||fje<p[-aV/(i- 


16. <C S C w.', v') 


a Sc u, V ) 
au 


^ = 2tt cr 


2 a W(cc, V) 

' a u 


S(LC-u') ^cv- v') 


17. 


< 


a 5 eu', v' ) a S ( V ) 


3 U' 


3U 


> = 2 'Tr <r. 




3LC^ 


S( M-u') Sc V- v'3 
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APPENDIX VII 


COX AND MUNK (1954) SUlsEACE MEAN SQUARE SLOPES 
1 . For clean sea 

+ CT* = 5.12 * lo + 0,003 ± 0,00.4 

cross ‘'P 

2 . For slick sea 

/T^ + <T^ = 1.56 * + 0.008 1 0.004 

‘^CI'OSS 'op 

In both (1) and (2), V is in meter/second. 

For isotropically rough sea surface, we hove 

'W. = .4. r O' + (j ) 

C,R a c- '.'O / 

where and are the rms slopes of the clean and the slick sea, respectively. 
Values of and ore given in Toble i for different svind speeds. 
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TAGLF I cox N 

MUVJK RMS SUREACE SL'O 

^ES 



N 

CLEAN SEA 

SLICK SEA __ 



v(m/se:) 

..... . Ml 

M2 

M3 


0.500 

0 . 2763E-01 
O,4o09E-(ii 

0.5273E-01 
0 . 6626E-ni 

0 . 5914E 
0 , 7 9 9 4 E - 0 1 

-01 

2.000 

0 . 6/97E-01 
0,5967E-o1 

0 . S136E-01 
0.7457E-01 

0.9204= 

0 . 3695E-01 

-01 

... '*.000 

0 .9H69P-01 
3.7155E-01 

. . .. 0.10P4E 00 

0 . 3439E-01 

0 . 1172E 
0 . 9550E-CI1 

0 0 

5.000 

0.1109= no 
0.76B1E-01 

0.1196E 00 
o.aaasE-oi 

0 . 1 2 7 7 F 
0.9950E-ni 

0 0 

6,0 00 . .. 

0 , 12t=3E 00 
r . . 3173E-C1 

0.129SE .3 0.,... 
n . 9317E-01 

... . 0 . 1373E 

0 . 10 3 3E nil 

L C 

8.030 

0.1414= 00 
0,90 7 7E-q,1 

0.14a3E 00 
.D.1012E CO- 

0.1549= 
. 0.11D6E 00 

0 0 

1 0,030 

0.15B4E-00 

0.9399£-ni 

.. .... .0 ,164 6E._0 0.. . 

0.10S6E 00 

, , 0 . 1 7 J 6 E 

0.1175E 00 

0 0 

12.030 

0.1738= 00 
3 , 1066.E .00. 

0.1795E 00 
..Q.1156E. 00 . 

0 .1650= 

. .,0.1239= DO . 

0 0 

14.0 00 

C . 1 8 -3 0 E 0 0 
0 . i 1 3 7 E no 

. Ot 3.9.3 2E 0 0., 
0.1221= 00 

0.1R53E 
0.1301= 00 

0 0 

15.000 

0 . i 9 4 7 £ no 
n . 1 1 7 0 E on 

0.1997E 00 
0.1253E CO 

n . 2 0 4 7 c 
0 . 1 .3 3 0 = 0 0 

n 0 

16.030 

n .2niiE no 
o'. 120 3= go' 

n . 206lE , 00. . 

0.12B4E 00 

0,2109= 
"O , 1 3 5 9 E 0 0 

0 

16.030 

0.2135= 00 
G.1266E OC 

n.PlBlE 00 
0 , 1343E no 

n . 2227G 
0,1416= 00 

0 0 

20.030 

0.2252E 00 

n.2296E 00 

0 .2339= 

DO 


G.1327E 00 

n,l400E 00 

0 . 1 4 7 0 E no 


23.000 

0.241 ftp 0 0 
0 . 1412E 00 

0.2457E 00 
0,l4BlE.00 

n .2496= 
. 0.1547E OC 

0 0 

25.030 

0.2520= 00 

n,2559E 30 

n . ?595E 

00 


0.1466E no 

0.1533E 00 

■ 0 . 1 5 9 7 E Q 0 

... i;....FOR EACH ',|riD Sp, The upper line is the HMS Sl'O 
N the lower line is roR SLi::K spa. 

2, M1,H2,N are TtiF 1. OWE R i ME D I JM , N D'^PER V.ALJES 

Sl“^eS, RespecTi vely. 

=>£5 =0R clean S 
OF RMS SURFACE. 

E A 1 
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